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Hereafter we shall present a brief overview on some of the most significant 
achievements obtained by means of NMR and /iSR techniques in highly frus- 
trated magnets. First the basic quantities measured by the two techniques will 
be presented and their connection to the microscopic static and dynamical spin 
susceptibility recalled. Then the main findings will be outlined, starting from 
the most simple frustrated units, the molecular nanomagnets, to artificially 
built frustrated systems as 3 He on graphite, to magnets with a macroscopically 
degenerate ground-state as the ones on a pyrochlore or kagome lattices. 



1 Some basic aspects of NMR and fiSH techniques 

NMR and /iSR are very powerful techniques which allow to investigate the 
microscopic properties of spin systems through the study of the time evolution 
of the nuclear magnetization M(t) and of the muon spin polarization P(t), 
respectively [HE]- Each technique has its advantages and disadvantages. In 
NMR one knows the crystallographic position of the nuclei under investigation 
and, therefore NMR results can be more suitably compared to theories. On the 
other hand, NMR experiments cannot be performed in compounds where just 
low sensitivity nuclei are present or where the fast nuclear relaxations prevent 
the observation of an NMR signal. Still, polarized muons can be injected 
into the sample and used as a probe of the local microscopic properties of 
the system under investigation. Moreover, by means of /iSR it is possible to 
detect relaxation times shorter than 0.1 [is, about two order of magnitudes 
shorter than the shortest relaxation time NMR can measure. Since the nuclear 
magnetization is the quantity detected in the NMR experiments, generally a 
magnetic field has to be applied to generate it. On the other hand, the muon 
beam is already polarized before entering the sample, so that the system under 
investigation can also be studied in zero field by means of /iSR. This aspect 
is particularly relevant if one wants to investigate the intrinsic properties of 
a certain system without perturbing it with a magnetic field. Nevertheless, 
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novel ground-states can be induced by the application of high magnetic fields 
(typically above 10 Tesla) where /iSR experiments cannot be performed while 
NMR experiment can. Hence, although both techniques appear to measure 
similar quantities (see next section), in view of the above considerations they 
are often complementary and their combination is a rather powerful method 
to investigate the local microscopic properties of frustrated magnets. 



1.1 Line shift and line width 

The time evolution of M(t) and P(t) is determined by the hyperfine interac- 
tions which can be summarized in the form 

H = H z + Hn-n + W„- c + HeFG ■ (1) 

The effect of all the four terms on the hyperfine levels and on the NMR 
spectra for 7 = 3/2 are depicted in Fig. [TJ We recall that the intensity of the 
whole spectrum is proportional to the nuclear magnetization. The first term 

Zeeman + n-n + n-e + quadrupole 
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Fig. 1. Schematic illustration of the modifications in the hyperfine levels of I = 
3/2 nuclei, due to the different terms of the nuclear hyperfine Hamiltonian. The 
corresponding modifications in the NMR spectra are reported at the bottom of the 
figure. 



describes the Zeeman interaction with an external field, while H. n -n is the 
dipole-dipole interaction among the nuclear spins or between the muon and 
the nuclear spins. This interaction yields a broadening of the NMR p] and /iSR 
[2] spectra. In certain compounds, as the cuprates, the nuclear dipole-dipole 
interaction is mediated by the electron spins and from the dipolar broadening 
information on the static electron spin susceptibility x'{Q) is obtained [3]. 
The last term TLefg describes the interaction between the nuclear electric 
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quadrupole moment Q and the electric field gradient (EFG) generated by the 
charge distribution around the nucleus. This term is non-zero when the nuclear 
spin I > 1/2 and is, of course, absent in the muon interaction Hamiltonian. 
The quadrupole interaction is very sensitive to the modifications in the local 
configuration and allows to evidence distortions induced by the spin-lattice 
coupling [4] or the presence of a non-homogeneous charge distribution induced 
by charge ordering, for instance [5]. The third term is the most relevant one 
to probe frustrated magnetism, as it describes the hyperfine interaction with 
the electron spins S. As most of the systems we shall be dealing with in the 
following are insulators, one can write 

Hn-e = -7Kj2liAikS k (2) 



with An- the hyperfine coupling tensor, I and S are the nuclear/yu + and elec- 
tronic spin operator respectively, i and k are the nuclear//Lt + and electron spin 
indexes, respectively, while 7 the nuclear gyromagnetic ratio. Then, the 
hyperfine field at the z-th nucleus//^ 4 " will be given by hi — J2k AikSk and in 
the presence of a non-zero average polarization < S > hi = J2 k Aik < Sk >■ 
Thus, one can directly estimate < S > from the precessional frequency 
uj = 7 J2k Ak < Sk > of the nuclei or of the muons around the local field. 
When an external field Hq \\ % is applied the local magnetic field becomes 

B = H a + Y,A k <S k > (3) 

k 

and the resonance frequency will be shifted to 

LU = LU L (1 + K) (4) 

with ujl = ^Hq the Larmor frequency and, for Hq ^> | J2k A k < S k > \, 



J2k A k <S k > 
K = ± '-± . (5) 

-no 

In general K is a tensor and K — J2k AkX{q = 0, w = 0). Hence, one notices 
that from the shift in the precessional frequency of the nuclei (or one can 
derive the static uniform susceptibility associated only with those electron 
spins which are coupled to the nuclei under investigation see FigO 

If, for some reason, the distance between nuclei or muon and the electronic 
spins is random or there are missing spins in the sample, the hyperfine coupling 
A will be a random variable leading to a distribution of precessional frequen- 
cies and to an increase in the relaxation of the polarization/magnetization. 
For an external field applied perpendicular to the initial polarization of the 
muons one finds 
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Fig. 2. Temperature dependence of 71 Ga(4f) NMR shift K in SrCrg p Gai2-gj,0ig 
(p = 0.95) compared to the macroscopic susceptibility derived with a magnetometer. 
Since Ga(4f ) nuclei are strongly coupled only to chromium ions within the kagome 
bilayer one can single out just the intrinsic susceptiblity of the latter, whereas the 
macroscopic susceptibility detects also the contribution from defects. [B] 



P x (t) =exp - 



cos(wi). 
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l/Tj could also represent the decay rate of the NMR signal after an RF pulse. 
Assuming a distribution of hyperfine fields in the z direction one can write 
Ak as a sum of a mean value Ak plus a fluctuating component 8Ak ■ For the 
distribution 

p(5A k ) = r— exp 



Mi 

2ol 



one finds that 



1/2 



is the width of the spectrum which has an average shift 

LU — 7-ffo 

= X _ 

k 



K = 



If (Tfe and A are temperature independent parameters we expect 

l/T* a: K 



(7) 



(8) 



(9) 



where the temperature is an implicit parameter. It is noticed that this pro- 
portionality is correct regardless of the form of the relaxation ([6]) . The break- 
down of the validity of Eq. [5] would indicate a modification in the hyperfine 
couplings, which is typically expected when lattice distortions take place (see 
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See's. 2.2 and 2.5). On the other hand, in certain systems, although the hyper- 
fine coupling is constant, the spin polarization can be site-dependent. Then, 
Eq[5] no longer holds and the line broadening and its shape reflects the distri- 
bution of the local spin polarization. |7J 

1.2 Nuclear and muon spin lattice relaxation rate 1/Ti 

The transitions among the hyperfine levels, driven by the time dependent 
part of the hyperfine Hamiltonian (not shown in Eq. Q}, modify the nuclear 
spin population on each level, namely the longitudinal component of nuclear 
magnetization. Thus, in a frustrated magnet from the time evolution of the 
nuclear magnetization it is possible to derive information on the spin dynam- 
ics, which drives the fluctuations of the hyperfine field. The recovery of the 
longitudinal component of the nuclear magnetization, after the nuclear spin 
ensemble has been brought out of equilibrium with an ad hoc RF pulse se- 
quence, is described by a characteristic decay rate 1 jT\ . In case of relaxation 
mechanisms driven by fluctuations of the hyperfine field h(t), by resorting to 
time-dependent perturbation theory (h is considered small with respect to Ho) 
and assuming that the frequency of the field fluctuations lu — 2ttv ^> 1/Ti, 
one can write 



This fundamental expression shows that 1/Ti is driven by the transverse com- 
ponents of the fluctuating field at the nucleus, owing to magnetic-dipole se- 
lection rules, and that 1/Ti is proportional to the Fourier transform of the 
correlation function at the resonance frequency, in order to allow for energy 
conservation. In other terms 1/Ti probes the spectral density at the resonance 
frequency ujl which is typically in the MHz range, orders of magnitude be- 
low the spectral range accessed by inelastic neutron scattering experiments. It 
should be noticed that this does not mean that from 1/Ti one cannot estimate 
relevant energy scales much larger than Tilul. In fact, when sum rules apply, as 
it is often the case in spin systems, the amplitude of the low-frequency spec- 
tral density is determined by the characteristic frequency of the fluctuations 

UJ » UJ L . 

In general, when collective spin excitations are present one can write 



and by substituting the transverse components of h(t) in Eq. 10 expression 
one finds that 




(10) 




(11) 




(12) 
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Fig. 3. Si form factor in the first Brillouin zone of the two-dimensional frustrated 
antiferromagnet L^VOSiO^ It is evident that excitations at wave-vectors (±7r/o, 0) 
or (0, ±7r/o) are filtered out, namely that 29 Si 1/Ti is not sensitive to these modes. 

One notices that, being the nuclei local probes, 1/Ti is related to the integral 
over the Brillouin zone of the component of the dynamical structure factor 
Saaiq^h) at the Larmor frequency. In Eq. 12 \A q \ 2 is the form factor which 
gives the hyperfine coupling of the nuclei with the spin excitations at wave- 
vector q (see Fig. [3]). The term _L indicates that one has to consider the 
products | A q \ 2 S aa (q, wo) associated with the perpendicular components of 
the hyperfine field at the nucleus. From the fluctuation-dissipation theorem, 
by recalling that usually fc^T 3> hu>L one can also write 



This rather general expression shows how the nuclear or muon spin-lattice 
relaxation is related to the spectrum of the excitations characteristic of each 
frustrated magnet. 

The above equations apply also to /xSR spin-lattice relaxation rate when 
a large magnetic field is applied. On the other hand, when /jSR operates in 
zero field, standard perturbation methods to analyze the data, such as Eq. llOl 
are no longer valid. This is because the transverse direction is not defined and 
the internal field is not small compared to H. Accordingly, different methods 
are required to account for the muon relaxation function in zero and small 
external fields. Moreover, usually one does not know the muon stopping site, 
and therefore the discussion is done using the field at the muon site B rather 
than the hyperfine coupling A. The treatment of the muon T\ in zero or small 
fields is done in two steps. The first step is the static case where T\ = oo. In 
the second step, the dynamic fluctuations are added and T± becomes finite. 




(13) 
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1.3 /xSR: the static case 

The fully polarized muon, after entering the sample, comes to rest in a 
magnetic environment. Since the mechanism which stops the muon is much 
stronger than any magnetic interaction, the muon maintains its polarization 
while losing its kinetic energy. Once the muon reaches its site, the muon spin 
starts to evolve in the local field B. The muon polarization P z along the z 
direction is given by the double projection expression 

P z (B,t) = cos 2 (9 + sin 2 (9cos(7jB|i) (14) 

where 9 is the angle between the initial muon spin and the local field direction 
(see Fig. HI). This angle is related to the field values by 

cos^ff , sin'*=^. 




Pz(0) 

Fig. 4. Muon spin polarization rotating around a magnetic field in an arbitrary 
direction. 



In a real sample, however, there will be a distribution of internal fields and 
the averaged polarization is 



P z {t) = / P (B) 



R2 P?2 , r>2 



d d B 



(15) 



where P z (t) is the sample averaged polarization, and p(B) is the field distri- 
bution which is normalized according to J p(B)dB 3 = 1. If the distribution 
of internal fields is only a function of |B| then we can write 

P z {t) = J p(\B\) [cos 2 6» + sin 2 6» cos( 7ai |B|<)] B 2 dBdQ. 



It is convenient to define p'(\B\) = 4wp(\B\), so that J p'(\B\)B 2 dB = 1 and 
the angular dependence can be integrated out giving 
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P*W = l + l J p'(|B|)cos( 7 JB|i)i? 2 di?. 

If, for example, the system has long range order the field at the muon site is 
centered around a specific value ujq/j^ so that 



2ttAB- 



■ exp 



B 



-I 



/2A 2 



then 



P z {u Q ,A,t) = - + -exp - 



A 2 f 



cos(ujQt) 



and oscillations will be observed in the data. At long time the polarization 
will relax to 1/3 since effectively 1/3 of the muons experience a field parallel 
to their initial spin direction and do not relax. 



H 




Pz(0) 

Fig. 5. Muon spin polarization rotating around the vector sum of an external mag- 
netic field in the initial muon polarization direction, and an arbitrary internal field. 



When a longitudinal field H is applied in the direction of the initial muon 
spin as in Fig. the situation becomes more complicated, and there is no 
closed form expression. However, some simplifications could be made to reduce 
the dimension of the integrals for the purpose of numerical calculations. For 
example, if the local field is completely random with a Gaussian distribution 
then 



P(B) 



7^ 



(27r) 3 / 2 /A 3 



exp 



2A 2 



In this case Eq. [15] could be simplified to [8] 



(16) 



P z (LU L ,A,t) = l- 



2A 2 



K) 2 

2A 4 rt 



1 - cxp(-i/\ 2 i 2 ) cos(w L i) 



K) 3 



f 1 

/ exp( — A 2 t 2 ) sin(wLT)dr 
Jo 2 



(17) 



This is known as the static-Gaussian-longitudinal-ficld Kubo-Toyabe (KT) 
function. Figure [6] shows P z (ujl, A, t) for a variety of ui^. Interestingly, despite 
the fact that the external field is in the muon spin direction, wiggles are seen 
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Fig. 6. Muon polarization function in a Gaussian internal field distribution and 
external field pointing in the initial muon spin direction. Different values of the 
external field H are shown. 

in the polarization, and their frequency is given by ui. When loi 3> A, the 
muon does not relax any more because the field at the muon site is nearly 
parallel to the initial muon spin direction. In this situation we say that the 
external field decoupled the muon spin from the internal field. Finally, in the 
zero field case (Hi, = 0) Eq. [17] reduces to [8] 

P z (0, A,t) = 1 + |(1 - AH 2 ) exp(~A 2 t 2 ). (18) 

This polarization function is known as the static-Gaussian-zero-field KT. At 
early time it has a Gaussian-like behavior. It reaches a minimum on a time 
scale set by A after which it recovers and saturates again to 1/3. 

1.4 /xSR: the dynamic case 

If we now add dynamics numerical methods must be applied. If the dynamic 
part of the local field at the muon site (5B fluctuates in time and magnitude 
so that 

3/\ 2 

(JB(i)JB(O)) = -=- exp(-2vt), (19) 
% 

where v is a fluctuation rate, and under the strong collision approximation, 
the muon polarization will obey the Volterra equation of the second kind 13J. 
The polarization P z (v, ujl,A, t), which now also depends on the fluctuation 
rate u, obeys 

P z (v,u) L ,A,t) = e- vt P~ z (0,u L) A,t) + 

vj dt'P z (v^ L ,A,t-i!)e- vt 'p z (b^ L ,Aj) (20) 
Jo 
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where P z (0, u>l, A, t) is the static relaxation function, namely, the polarization 
if the local field was frozen in time. The factor e~ vt is the probability to have 
no field changes up to time t. The factor e~ ut vdt' is the probability density to 
experience a field change only between t' and t'+dt'. The first term on the r.h.s 
is the polarization at time t due to muon that did not experienced any field 
changes. The second term on the r.h.s is the contribution from those muon that 
experience their first field change at time t' . The factor e~ ut P z (0, H, A, t')vdt' 
is the amplitudes for the polarization function that evolves from time t' to t, 
which can involve more field changes recursively. This equation can be solved 
numerically |10j . 




2 4 6 8 10 

At 

Fig. 7. Expected muon relaxation in a dynamic field with Gaussian instantaneous 
distribution and no external field. Different values of fluctuation rates are shown. 

There are three ways of using Eq. [201 to obtain dynamic information. The 
first is in simple cases where -P z (0, u>l, A, t) is known analytically as was done 
by Brewer et al. [TT] for F-/1-F bond. The second is when P z (0, u>l, A, t) 
must be obtained numerically as in the cases of Gaussian [8] or Lorenzian 
[12] field distribution with external longitudinal field. The third way is to 
measure P z {0, vl, A, t) by cooling the system to low enough temperatures 
that dynamic fluctuations are no longer present, and to use the measured 
P z (0, ujl, A, t) as an input to the Volterra equation [T3"] . 

Taking the polarizations generated by static field distribution given by 
Eq. [17] with lul = 0, and using it as an input in the Volterra equation, gives 
the dynamic polarizations shown in Fig. [7] We see from this figure that when 
dynamic fluctuations are present the 1/3 recovery is lost. As the fluctuation 
rate v increases, the Gaussian-like relaxation at t — > is also lost, and when 
v > A the relaxation becomes exponential. 
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Finally, in Fig. [5] we present the most complicated relaxation function 
combining Gaussian field distribution, fluctuations, and longitudinal field. 
This is known as the dynamic-Gaussian-longitudinal-field-KT relaxation func- 
tion. We have chosen a special value of the parameters A = 11.8 Mhz and 
v = 12.2 MHz for reasons that will become clear in Sec. 2.5, and show how 
the polarization is modified as H varies. 




0.0 0.1 



0.2 0.3 0.4 
Time (\i sec) 



Fig. 8. Solid lines: Expected muon relaxation in a combination of internal field fluc- 
tuations with Gaussian instantaneous distribution, and longitudinal external field 
H. Symbols: Relaxation data from SrCr9 P Gai2-9 P Oi9 from Ref. [42] in two different 
values of H. The disagreement between data and model indicates unusual behavior 
and is discussed in Sec. 2.5 



It is interesting to mention that for the case v > A there is an approx- 
imate expression for the dynamic-Gaussian-longitudinal-field-KT relaxation 
function [14j given by 

P z (t) = exp(-r(t)t) (21) 

where 

2 A 2 {[lo 2 + v 2 ]vt+[uo 2 - v 2 ][l - e~ vt cos(w L <)] - 2vu^e~ vt sm(w L i)} 



r(t)t 



K + ^ 2 ) 2 



In the long time limit (vt ^> 1) one finds that lim^oo P z (t) — Pq exp (— t/Ti) 
leading to the standard expression for Ti 

1 2A 2 v 



Ti K + z/ 2 ) 



( 22 ) 
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This expression demonstrates that when the external field is small, namely, 
uj l <C v, T\ will show no field dependence. But when the field is large ui L ^> u, 
T\ will increases with increasing field. Therefore, field dependent measure- 
ments can be used to provide information on v, and distinguish between static 
and dynamic cases. 

Equation l22l could be related to Eq. [T3l given for NMR \/T\. In systems 
characterized by a spin-spin correlation function decaying in time as exp(—vt), 
neglecting the q— dependence, one can write for the spin susceptibility 

nv kb-1 v + uj z 

so that A 2 in Eq. 22 can be related to the amplitude of the spin fluctuations 

4 2 = gs(s + i)£(KiV 

When the applied field is small or the fluctuations are fast Eq. [22] reduces to 

1 ^ (23) 



Tj v 

Typical values are T\ ~ 0.1 /xsec, A ~ 10 MHz, and v ~ 10 /zsec -1 . 



2 From zero to three-dimensional frustrated magnets 
2.1 Molecular Magnets 

In recent years major attention has been addressed to the study of molecular 
crystals containing molecules formed by paramagnetic ions with significant 
intramolecular exchange couplings and negligible intermolecular couplings, so 
that each one of them can be considered as an independent nanomagnet [15] . 

In some nanomagnets V 4+ S — 1/2 ions form a triangular lattice. This 
is the case, for instance, of V15, where at low temperature (T) 12 V 4+ spins 
are coupled in singlets and the remaining 3 ions form a triangle. From the 
study of 51 V NMR shift below 100 mK Furukawa et al.[T|)] have preliminarily 
estimated the hyperfine coupling and then estimated the expectation values 
for V 4+ magnetic moments in zero-field. The resonance frequency of the i-th 
nucleus is given by Wj = "j(H +A < Si >), where the internal fields A < Si > 
have a different orientation and magnitude, depending on < Si >. In order to 
derive < Si > the authors have studied the magnetic field dependence of a;,, 
as shown in Fig. [5J The expectation values were found to be consistent with 
a doubly degenerate ground-state of the form 
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Fig. 9. Magnetic field dependence of the resonance frequency for different 51 V 
sites in Vig, used to derive the spin polarization of each V 4+ |16| . In the inset the 
schematic view of the twofold degenerate ground-state is shown. 

Ipa = -j=(\ HT> -| ITI>) 

= ^(| TU> -I m> -I UT>) (24) 

where f or J, represents the orientation of each one of the three spins. /iSR lon- 
gitudinal relaxation rate measurements in the same compound [17] evidenced 
a nearly T-independent relaxation rate at low temperature, which is possibly 
associated with the transitions among these two degenerate ground-states. If 
one considers that the frequency of the spin fluctuations lu 3> lol > one can de- 
rive the characteristic tunneling rate between the two degenerate states. We 
point out that a nearly T-independent relaxation is observed also in frustrated 
magnets with a macroscopically degenerate ground-state (see Sect. 2.6). 



2.2 Antiferromagnets on a square-lattice with competing 
interactions: the Ji — J2 model 

V 4+ ions can also form other magnetic structures characterized by a strong 
frustration of the magnetic moments. In fact, certain vanadates can be consid- 
ered as prototypes of frustrated magnets on a square lattice where frustration 
arises from the competition between nearest neighbour (Ji) and next-nearest 
neighbour {J 2) exchange couplings. The first NMR studies have been carried 
out in Li 2 VOSi04. 7 Li NMR spectra were observed to split in three different 
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Fig. 10. (circles) Temperature dependence of the order parameter in the collinear 
phase of L^VOSiO/i, derived from zero-field /xSR. (squares) Temperature depen- 
dence of 7 Li NMR linewidth for H \\ ab in Li2VOSi04. It is noticed that a broaden- 
ing starts well above the transition temperature (vertical dotted line) and is possibly 
due to a frustration driven lattice distortion. 



lines for T < T c ~ 2.9 K [T51 one unshifted and the other two sym- 
metrically shifted with respect to the central one. This splitting of the NMR 
line was the first evidence that this compound is characterized by a magnetic 
collinear ground-state, as confirmed few years later by neutron scattering ex- 
periments [20] ■ A careful study of the order parameter has been carried out 
by means of zero-field /iSR measurements (Fig. [T0|). where the /i + polariza- 
tion is characterized by oscillations at a frequency directly proportional to 
V 4+ magnetic moment [5T] (Sect. 1). The continuous increase of the order 
parameter for T — ► T c was found to be described by a critical exponent close 
to the one expected for 2D XY universality class. Above T c , when no internal 
field is present, the NMR shift K is expected to be proportional to the static 
uniform spin susceptibility Xs (see Eq. 5), with a slope given by the hyperfine 
coupling. Remarkably, for T c < T < J\ + J 2 a change in the slope of 7 Li K vs 
Xs is noticed [22] , suggesting a modification in the average hyperfine coupling 
(see Sect. I). Moreover, the broadening of 7 Li NMR linewidth in the same 
temperature range can be explained only if the distribution of the hyperfine 
couplings is also increasing in the same temperature range. We remark that in 
this temperature range the average shift decreases on cooling below 5 K while 
the linewidth increases, namely the breakdown of Eq. 9 is noticed. A similar 
scenario is observed for 95 Mo NMR in M0VO5 [23] • These modifications in 
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the hyperfine coupling have been associated with a lattice distortion driven 
by the spin-lattice coupling, which releases the degeneracy of the magnetic 
ground-state. 
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11. Comparison of the temperature dependence of 7 Li 1/Ti in Li2VOSi04 
63 Cu 1/Ti in La2Cu04, where frustration is negligible. The temperature is 



normalized to Curie- Weiss temperature 0, while 1/Ti is normalized to its value at 
T»6>. 



The T-dependence of the in-plane correlation length £ of Li2VOSi04 has 
been estimated from 7 Li 1/Ti [Mj- In fact, by resorting to scaling arguments 
one can rewrite the dynamical structure factor in Eq. 12 in terms of powers 
of £. One finds 

^ ~ T ~ exp(2^ Ps /T). (25) 
J-i 

with p a the spin stiffness and z the dynamical scaling exponent which is 
estimated around z — 1 on the basis of some physical considerations [21] . 
It turns out that for a two-dimensional antiferromagnet on a square-lattice 
frustration yields a less pronounced increase of £ on cooling, namely to a 
decrease in the spin stiffness. 

The decay of the longitudinal muon polarization in Li2VOSi04 and 
Li2VOGeC>4 was observed to be well described by Eq. 18 and evidenced a 
dynamics at frequencies well below the Heisenberg exchange frequency [24] . 
This dynamic has been ascribed to the fluctuations within the two-fold de- 
generate ground-state over a barrier [25j 

em . 
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2.3 Magnetic frustration on a triangular lattice 

The simplest two-dimensional lattice frustrated by the geometry of the in- 
teractions is the triangular one. Several systems on a triangular lattice have 
been investigated in the last decade, either insulating [26] or metallic ones 
|27] . Some of these compounds display a rather interesting phase diagram as 
a function of the magnetic field intensity, as will be shown in a subsequent 
chapter of this book by M. Takigawa and F. Mila. Also in re-(ET) 2 Cu 2 (CN)3 
molecular crystals, which are insulators at ambient pressure, the spins are 
arranged on a triangular lattice. Recently, it has been observed (Fig fT2")) that 
1 H 1/T\ abruptly decreases at low temperature [28 . This decrease can be 
associated with the onset of a gap A between a collective singlet and triplet 
states, which yields 1/Ti oc exp(— A/T) for T A, and can be considered 
the first evidence for a spin-liquid phase in a triangular antiferromagnet. This 
could be possibly the experimental evidence of the long sought RVB state 
predicted by Anderson long ago [29] . Interestingly enough the application of 
hydrostatic pressure was observed to drive this system into a superconducting 
ground-state which can still be justified in terms of an RVB description. 
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Fig. 12. 1 H nuclear spin-lattice relaxation rate 1/Ti for a single crystal (open 
circles) and a poly crystalline sample (closed circles) of k-(ET)2Cu2(CN)3 and a 
single crystal of jt-(ET) 2 Cu[N(CN) 2 ]Cl (open diamonds) [28]. The Inset shows the 
low temperature part of the data in logarithmic scales. 
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Another spin system on a triangular lattice which has attracted much 
attention in recent years is Na^CoC^, where a rich phase diagram develops 
upon Na doping [3U]- An accurate study of the 23 Na and 59 Co NMR spectra in 
oriented powders evidenced a charge order for x ~ 0.7 [5]. The signature of the 
charge order is the presence of three distinct Na sites characterized by different 
quadrupole couplings and magnetic shifts, which imply a well defined order of 
the Na + ions and of the Co charges in the C0O2 planes. On the other hand, 
for x = 0.5 the study of the temperature dependence of 59 Co NMR spectra 
reveals that the the electric field gradient at the Co site does not change at the 
metal-insulator transition, indicating the absence of any charge ordering 31j. 
These NMR measurements have allowed to clarify the nature of the ground- 
state of this system in the doping range 0.5 < x < 0.75 [32] • When full Na 
doping is achieved the system eventually becomes non-magnetic and the NMR 
shift vanishes [33]. Remarkably for x — the isostructural compound C0O2 
is found to be metallic, with a crossover from a strongly correlated metal to 
a Fermi-liquid behaviour at low temperature [34j . In fact, around 4 K the 
temperature dependence of \jT\T starts to flatten, as expected in a Fermi- 
liquid pp. 




T(K) 

Fig. 13. From the top: specific heat, muon spin-lattice relaxation rate and 23 Na 
NMR signal intensity in NaCrCh vs. temperature. It is noticed that the peak in the 
specific heat and 1/Ti occur at different temperatures. The loss of 23 Na signal is 
possibly associated with an increase of the NMR relaxation rates. 
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Recently, it has been pointed out that NaCr02 is an excellent realization 
of a S = 3/2 triangular Heisenberg antiferromagnet. Remarkably, in this com- 
pound while specific heat and magnetization measurements indicate the onset 
of a transition around T c ~ 40 K, both muon spin rotation and NMR reveal 
a fluctuating crossover regime extending well below T c , with a peak in 1/Ti 
around 25 K [35] (see Fig [T3"|) . This apparent discrepancy might indicate the 
presence of vortex-antivortex excitations decoupling around 25 K. 

Magnetic frustration can be associated not only with the geometry of the 
electron spin arrangement but also with the one of indistinguishable nuclear 
spins, as it is the case for 3 He. The triangular lattice topology can be achieved 
by evaporating a single layer of 3 He on a graphite substrate. Then one can 
conveniently use the intensity of the NMR signal, which is proportional to 
the nuclear magnetization, to track the T dependence of the nuclear spin 
susceptibility : 36J. 3 He NMR measurements have been carried out down to 
tens of fiK, more than an order of magnitude below the exchange coupling. 
The low temperature increase in the nuclear spin susceptibility evidenced the 
relevance of higher order multiple spin exchange interactions which, together 
with the triangular geometry of the nuclear spins, makes the system strongly 
frustrated and possibly characterized by a gapless spin-liquid ground-state 

|33 ESI- 
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Fig. 14. The muon shift K and transverse relaxation time I/T2, versus temperature. 
Inset, a plot of I/T2 versus the system's suspetibility. 



2.4 /xSR and NMR in the spin-1/2 kagome lattice ZnCu 3 (OH) 6 Cl 2 

A very interesting model compound for the study of magnetism on the kagome 
lattice is the herbertsmithite ZnCu3(OH)gCl2. The moments in this material 
originate from Cu 2+ which has a spin 1/2. Therefore, it is ideal for the in- 
vestigation of quantum ground states. Unfortunately, different probes such as 
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muon [3H], O, [35] Cu, and CI [40j nuclear magnetic resonance suggest dif- 
ferent behavior of the shift below ^ 50 K and the origin of these variations 
are not clear yet. Here we present only the /xSR results. First we examine the 
muon shift and T 2 * in ZnCu3(OH)6Cl2 independently (Fig. [T4"]) . and also one 
as a function of the other (inset of Fig. [T4]) . The shift, and hence the sus- 
ceptibility, increase continuously upon cooling and saturate below 200 mK. 
This indicates that the spin- 1/2 kagome does not freeze or form singlets. The 
ground state is paramagnetic. In addition, the K and 1/T 2 * follow each other 
as the temperature is lowered, as shown in the inset, as expected from Eq. 
[9] Although K is not exactly a linear function of 1 /T 2 * there is no reason to 
suspect a modification in the hypcrfine coupling upon cooling, namely that 
lattice deformation is present in this case. 

Second we examine whether the ground state is separated by a gap from 
the excited ones. If such a gap exists it will take a finite temperature to 
generate excitations and achieve a non zero x" QQ (q, ujq). Therefore, according 
to Eq. [T3l 1/(T]T) should extrapolate to zero. This has not been observed 
experimentally gU gQl [39] . The 37 CI 1/Ti depicted in Fig. [15] is proportional 
to T meaning that x" aa(q,wo) is finite in the T — > limit. Therefore, there 
is no evidence of a gap and the kagome lattice seems to be an exotic magnet 
with no broken continuous symmetry but gapless excitations. 




2.5 The problem of the /x+ relaxation in some kagome lattices 

At the very beginning of the research in the field of frustrated magnets it was 
noticed that the muon relaxation function is unusual [42]. The symbols in Fig. 
[8] show the polarization at a temperature of 100 mK in the kagome system 
SrCrg p Gai2_9pOi9 (SCGO) in zero field and a longitudinal field of 2 kG [42] . 
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First, no oscillations are found, so the internal field is random with either 

static or dynamic nature. Second, the relaxation at early time is Gaussian, 

with a time scale of 0.1 /usee, so A must be on the order of 10 MHz. Third, 

there is no recovery so there must be some dynamic as in Fig. [7J But it must 

be that v ~ A. Had v been larger then A, the initial relaxation would have 

been Lorenzian (See Eq. I2T1 to l22l and Fig. [7])- Had v been much slower, the 

polarization would have recovered at least partially. In these circumstances a 

field of 2 kG which is equivalent to ojl = 170 MHz should have "decoupled" 

the relaxation. This is not happening. The solid lines in Fig. [5] represent the 

expected decoupling which is very different from the observed one. A model 

has been proposed to explain this problem [42 , which received the name 

sporadic dynamic (SD) [43] [44]. 

In this model, the sample is not relaxing the muon spin all the time 

but only for a fraction / of the time, as demonstrated in Fig. [111 In zero 

field it is clear that such a case will lead to sporadic dynamic polarization 
s d 

P z (v,0,A,t) = P z {v, 0, A, ft). However, even when the H is applied, the 

polarization changes only when the internal field relaxes the muon spin. In 

other words, the flat parts in Fig. \W\ stay flat even when H is on. Therefore 

we expect P| d (^, cjl, A, t) = P z {v, wl, A, ft) for all values of wl- Since u, cjl, 

and A always enter the relaxation function as a product with t (for example 

sc i 

see Eq. [5T]) we must have P z {v, wl, A, t) = P z {fv, fuj L , fA, t). 




o 2 4 6 8 10 

Time 

Fig. 16. Muon relaxation function in a case where the muons relax only sporadically 
at certain intervals of time. 



When analyzing the data we are actually estimating fA and fv to be 
10 MHz and = P 2 (10, fuj L , 10, t). Therefore, the effect of the H is reduced 
by factor /. In the fraction of the time the field is turned on, A and v are 
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much higher than 10 MHz. This is the reason we do not see decoupling. This 
model is very successful in explaining fiSK data. 




0.1 1 10 100 

Temperature (K) 



Fig. 17. Temperature dependence of the muon relaxation rate 1/Ti at H=50G. In- 
set, the rate at which Tj varies with field as a function of magnetic ion concentration 
p normalized by the percolation threshold p c . 

The model could lead to two interpretations. The first one is that most 
of the time the field at the muon site is zero due to singlet formation but 
every now and then a singlet breaks giving rise to a fluctuating field and 
relaxation 42J. The problem with this interpretation is that the muon relax- 
ation is temperature-independent below 2-3K, similar to the case presented 
in Fig. 1171 This usually happens when the system is in its ground state. But 
in the ground state there cannot be time evolution, namely, the field cannot 
turn on and off and the muon spin cannot relax. If the system is not in the 
ground state at T ~ 2 K, then a lower energy scale, below about 100 mK, 
must exist which separates the ground from the first excited states. Accord- 
ingly, the constant relaxation observed for lOOmK < T < 2 K, would arise 
from the quantum fluctuations within a manifold of weakly coupled nearly 
degenerate states. Another possibility is that the muon is hopping, although 
not evenly, between two sites with different relaxation rates. This, however, 
is also unusual at such low temperatures. A full description of the muon spin 
behavior in SCGO and similar systems is still lacking. 
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2.6 Persistent dynamics and lattice distortions in the pyrochlore 
lattice 

If there is one common message deducible from the study of frustrated mag- 
nets by /iSR it is that these systems maintain a fluctuating part of the moment 
even when the temperature is lowered much below the coupling energy scale. 
This is manifested in the saturation of 1/Tj upon cooling. In zero field 1/Tj 
is given by Eq. [53] Therefore, for 1/Xi to remain finite, at least part of the 
moment must continue to fluctuate and A cannot be zero. While both T\ and 
A can vary from system to system, the typical scale for v is /isec -1 . 

In Fig. [T7] we present what might be the canonical example of persisting 
fluctuations to the lowest temperatures, Tb 2 Ti 2 C>7 [46] . The full black symbols 
denote l/Xi over a wide temperature range. The increase in 1/Ti upon cooling 
indicates that v decreases, namely, the spin fluctuations slow down. However, 
in standard magnets, at some temperature long range order or spin freezing 
sets in, and the amplitude of the fluctuations 5B (see Eq. [19]) also decrease 
upon cooling. This is manifested in a decrease in A 2 upon cooling. The net 
result is a peak in 1/Ti at, or close to, the critical temperature. The 1/Ti peak 
is missing in Tb 2 Ti 2 C>7 suggesting that no static magnetic moment develops 
in this system. Similar results where obtained in other pyrochlore lattices such 
as Tb 2 Sn 2 7 07] and Gd 2 Ti 2 7 @4]. 

This conclusion leads to yet another open question in this area of research, 
namely, which type of excitation will dominate at low temperatures: spin wave, 
spinless, or spinons? This can be examined by measuring 1/Ti as a function 
of magnetic ion concentration p above and below the percolation threshold 
p c . A strong dependence of 1/Ti on p close to p c would suggest that the 
fluctuations emerge from a collective phenomenon. In contrast, if 1/Ti varies 
smoothly across p c it would suggest that the excitations are local in nature 
and impartial to the coverage of the lattice. In Fig. [T7] 1/Ti is presented for 
(Tb p Yi_p) 2 Ti 2 07 samples in which the Tb magnetic ion is replaced by non- 
magnetic Y [49]. Clearly the fluctuations have similar behavior both above 
and below the percolation threshold which is at p c — 0.39 [45] , 

Moreover, the muon Ti was found to be a linear function of H [49] . This 
is a very different behavior from Eq. [2H and suggests that the field correlation 
does not decay exponentially but rather with a power law. In this case, the 
dynamic properties of the system could be collected into dTi/dH . In the inset 
of Fig.[17]we show dT\/dH as a function of p. This quantity shows no anomaly 
at Pc, suggesting that Ti is controlled be local excitations. 

Finally, as in the case of Li 2 VOSi04 (see Sec. l2.2[) . magnetoelastic coupling 
is a very important interaction [50] [5T] in the pyrochlore system as well. It 
could lead to a lattice deformation in order to relieve the frustration. In fact, 
regardless of how small this interaction is, it will cause some lattice deforma- 
tion. In other words, there is no such thing as a perfect Heisenberg pyrochlore 
lattice since it must distort. Of course, for very small magnetoelastic coupling 
the distortion could be undetectably small. Nevertheless, the investigation 
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of the ground state in the presence of magnetoelastic coupling is a growing 
theoretical sub- field accompanied by experimental search for this effect. 

As explained above, a comparison between the muon relaxation rate and 
the muon shift or macroscopic susceptibility can provide a hint about lattice 
deformations. For example, in the pyrochlore Y2M02O7, 1/T| depends expo- 
nentially on the susceptibility. This is demonstrated in Fig. [TH] on a semi-log 
scale. The fact that the lattice parameters in Y2M02O7 vary upon cooling 
received confirmation also from NMR |51j and x-ray [52] experiments. In con- 
trast, the lattice of Tb2Ti2C>7 does not distort [53], and, indeed, 1/T 2 * depends 
linearly on the susceptibility as shown in the inset of Fig. [TBI 
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